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Abstract 

A simple queueing approach for segregation of agents in modified 
one dimensional Schelling segregation model is presented. The goal is 
to arrive at simple formula for the number of unhappy agents remaining 
after the segregation. 

The importance of the original Schelling model of segregation was mainly 
due to the combination of simplicity and ability to reproduce phenomena 
observed in sociology. Particular success of the model was related to descrip- 
tion of residential segregation. Whether for a street or a city, the simple 
ID or 2D (checkerboard) geometries were reasonable representations of real 
situation. And the results have proven to enhance our understanding of the 
social behaviour. 

It is worth noting that since 1990's considerable amount of literature 
pointed out that there are many aspects of social phenomena which may be 
described in dynamical agent models, but which require going beyond simple 
nearest neighbour ID and 2D models. Several phenomena, for example clique 
formation and clustering are much better represented in the network based 
approach. Network models allow also much greater flexibility in describing 
interactions between agents, in particular showing fat-tailed distribution of 
the number of neighbours and 'small world' phenomena. In general, models 
based on nearest neighbour ID and 2D geometries can provide only limited 
range of description of real life phenomena, where interactions between agents 
are usually long range and structureless. 

However, even such limited models may provide a glimpse of understand- 
ing of the mechanisms and processes. With this goal in mind, this paper 
presents a very simple approach to dynamics of segregation of agents in ID, 
short interaction case. 



The model is based on assumptions similar to those published in |T]. We 
consider a one dimensional chain of L sites and N agents (L > N, both 
being large), who can be in two states, + and — . Some of the sites in a chain 
are vacant. The densities of the +, — and vacant sites are considered to be 
constant and denoted, respectively, by p + , p^ and po- To stress the similarity 
of the model to spin dynamics it is useful to introduce another quantity, 
corresponding to overall magnetization density, defined as m = p + — p_. 

An agent is considered to be unhappy when certain fraction /* of its 
neighbours is of different type. We can model in this fashion the reactions 
of agents to 'unfriendly' neighbourhood. In the simplest model, considered 
here, we look only to the nearest neighbours as influencing the agent's be- 
haviour. This behaviour is simple: an agent may try to improve its situation 
by moving to a vacancy that would improve (or at least not worsen) its sit- 
uation. Starting from a random distribution of the agents the subsequent 
jumps should lead to Schelling segregation of the population. 

We shall further limit ourselves to the simple case of the 'constrained 
model' where only unhappy agents are motivated to move. The 'uncon- 
strained' situation, when also happy agents move as long as this does not 
make them unhappy, introduced by [1] is not considered here. The basic 
quantity we are interested in here is the number of unhappy agents at the 
end of the process of jumps, Uoo. Instead of random selection of the agents 
considering whether to jump or not we shall use a simple queueing model, 
in which agents play a sort of 'musical chairs' game. We assume an 'intel- 
ligent', God's-eye view algorithm designed to maximize the results of the 
optimization process. Such model of segregation presented here is rather 
close to the original Schelling algorithm, which has assumed that the un- 
happy agents would be queued up and attempt to improve their situation in 
ordered fashion. 

In the following we will use, as far as possible, the notation of [TJ. Let's 
assume without decreasing the generality of the treatment that m > 0, i.e. 
p+ > P-. We would thus call + agents the majority agents, and the — would 
be the minority. 

For the ID nearest neighbour model a given agent may have 0, 1 or 2 
agents as neighbours. As noted before, /* = 1/2 condition means that an 
agent is unhappy only if it is surrounded by two opposite agents. Let us 
consider this value of /* first. 

For a random initial configuration of agents we can calculate probabilities 
of + and — agents to be in unhappy state. 
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This can be expressed through m and po : 

P+ = ((l-po) 2 -m 2 )(l-p -m)/8 (3) 
p u _ = ((1 - p ) 2 - m 2 ) (1 - + m) /8 (4) 

It is worth noting that the average initial number of unhappy minority (— ) 
agents is greater than the number of unhappy majority agents. 

Starting from such a random configuration let us now try to estimate 
the outcome of the dynamic process, in which unhappy agents improve their 
utility function. In the constrained model only unhappy agents are allowed 
to move, and only if the move improves their utility function. 

Let's describe here the queueing approach. Let us assume that somewhere 
in the system there is an initial vacancy 'friendly' for a + agent. All that 
is needed at this stage is a single such vacancy. (It is worth noting that 
in the original Schelling model, and agent was simply inserted between two 
other agents, regardless of any formal vacancies). We move an unhappy + 
to this vacancy. But for a + to be unhappy it must have come from a — I — 
environment. And thus it leaves behind a vacancy (—0—), which is for sure 

'friendly' for an unhappy — agent (from a H h environment). The next 

step is then migration of such unhappy — agent, which, in turn leaves a + 
friendly vacancy (+0+). The process would be repeated until we run out 
of the suitable candidates for swapping, i.e. it would swap (into happiness 
state) p+ + and — agents. It is worth noting that this process is feasible 
even if there is just one suitable initial vacancy. In such case, even a random 
process would unerringly pick the right 'targets'. 

The remaining unhappy minority agents, p v _ — p+ could improve their 
situation if there are suitable vacancies in the system. The number of such 
'— ' friendly vacancies is 

p F v - = p (l - pi) = p Q (l - (1 - p + m) 2 /4) (5) 

The final amount of unhappy agents (all of minority type) would thus be, 
on average, 

Moo = P-~ P+- Py~ ( 6 ) 

Figure 1 shows the behaviour of for some values of m as function of the 
vacancy density po. It is worth noting that with the queueing method, if m = 
and po — ► the amount of unhappy agents for a specific starting condition 
would be proportional to statistical difference between initial distribution of 
unhappy + and — agents. This would generally be proportional to y/Z, so 
in the limit of large L should become negligible. The average value for 
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many initial configurations for m = is zero, as it is in the Dall'Asta et. al. 
unconstrained model. 

In the /* = 1/2 case we can compare the results of the queueing model 
with Figure 1 of the Dall'Asta et. al. paper p. The results are qualitatively 
similar, for example they show the same gradual decrease of the density of un- 
happy sites with increasing p and 'perfect happiness' above some threshold 
value Po( m )- The absolute values are, however bigger roughly by a factor of 2, 
which reflects the simpleness of the current model (for example the fact that 
one transition may make more than one agent happy, for example moving 
out the central + agent from the environment (H 1 h) makes, simul- 
taneously, both neighbouring — agents happy. This should further decrease 
the amount of unhappy agents u^). Such considerations were neglected in 
present reasoning for the sake of simplicity. 

We now turn to situation of /* < 1/2 when the process changes pro- 
foundly. Now an agent is happy only if it is surrounded by the same type 
agents. Thus, in the constrained model migration is possible only to vacan- 
cies that do not have any opposite agent in the neighbourhood. For example, 
to migrate a + one needs a vacancy of three following types: +0+ or 00+ or 
+00. All other vacancies are unsuitable. The ability to successfully migrate 
is thus much lower. 

When we start, as before, from a random ID chain we can calculate the 
initial number of unhappy agents 



It is much more difficult to build a successful queueing system now, as 
it is no longer true that a successfully migrating unhappy agent would leave 
behind an environment ready for accepting the opposite agent. In fact only 
in — h 0, OH — and — I — would migration of the + create a vacancy ready 
to accept a — . Let's call such positioning 'optimal'. 

The queueing procedure would then be divided into three steps. First, 
we would use a single 'starting vacancy' to make the chain of exchanges, but 
only between the optimally placed agents. The densities of such optimally 
placed agents is different for + and — and is given by 




(7) 
(8) 



or 
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P- = P+(2poP-+p') 



(12) 



where s stands for 'swap friendly' density, for example p s + is the density of 
— sites that would leave a friendly vacancy for a +. It is worth noting that 
p s + > p s _. As the exchanges go one by one, only p s _ of them can be made. 
This step reduces the number of unhappy sites to 

Pi - P U + -P- (13) 
P U - - P U -P- (14) 

The next step is to look for vacancies already present in the system that 
could accommodate, in a single transfer, an unhappy + or — . Such vacancies 
are, respectively, given by +0+, +00, 00+ and —0—, —00, 00— configurations. 
Their average densities are given by 

Pv = Po(2p p+ + P+) (15) 
Pv = Po(2p P-+p 2 ) (16) 

The remaining (if any) unhappy agents of any type can then, at the last 
stage, be moved to totally empty environments (000), which have average 
density of pvvv = Po- Thus the final density of unhappy sites is 

Moo = p+ + p U - - 2pi - pv - p v ~ Pvvv (17) 

Figure 3 shows examples of Uoo for some values of m. In the /* < 1/2 
model, the average number of unhappy agents at po is 1/2. Figure 4 
shows dependence of the on m for a small value of p . 
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r = 1/2 




Figure 1: as function of p for several values of m for the /* = 1/2 case. 
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r = 1/2 
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Figure 2: Uoo as function of m for small value of p = 0.001 for the /* = 1/2 
case. 
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r < 1/2 




Figure 3: as function of p for several values of m for the /* < 1/2 case. 
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r < 1/2 




Figure 4: Uoo as function of m for small value of p = 0.001 for the /* < 
case. 
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